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We review two commonly used approximate treatments of beamstrahlung at lepton colliders. We discuss the appli- 
cability of those approximations and show that they are bound to fail at very high energy e"*"e~ colliders such as 
CLIC. We model the beamstrahlung effects at CLIC by an exact numerical solution to the evolution equation for 
beamstrahlung and compare to the two approximations. We also discuss the impact of beamstrahlung on precision 
measurements of new physics parameters. As an illustration we consider Kaluza-Klcin lepton production in Universal 
Extra Dimensions and study the deterioration of the endpoint lepton energy spectrum caused by beamstrahlung. 



1. INTRODUCTION 

The fixed center of mass energy Ecm is a ceiebrated virtue of lepton colliders. In contrast to hadron colliders, 
where the parton- level Ecm varies from event to event, the additional kinematic constraint allows for a number of 
new measurements. Some well-known examples are: the presence of sharp endpoints in the distributions of various 
kinematic observables, the possibility of threshold scans, measurements of the missing mass, the angular distributions 
at production, etc. The availability of all these new experimental tools allows for precise determination of new physics 
parameters and is the primary motivation for pursuing the construction of a next generation lepton collider. 

One should keep in mind, however, that in reality a fraction of the beam energy is always lost due to the effects 
of beamstrahlung and bremsstrahlung, so that the true center of mass energy in the event is often displaced from 
its nominal value. The purpose of the present note is to investigate to what extent these effects can represent a 
problem for precision measurements of new physics parameters at very high energy e'^e~ colliders such as CLIC. To 
this end, in Sec. |21 we first review the two commonly used parametrizations of beamstrahlung. Then in Sec. O we 
discuss the applicability of those approximations for lepton colliders with sub-TeV Ecm (such as ILC) as well as 
Ecm in the multi-TeV range (such as CLIC). We find that the latter case falls outside the range of validity of the 
two approximations, and an exact treatment of beamstrahlung effects is required. We solve numerically the evolution 
equation for beamstrahlung using the exact expression for the Sokolov-Ternov transition probability (see Sec.|21l. In 
Sec. 0] we consider the example of Kaluza-Klein muon production in Universal Extra Dimensions (UED) and use 
our numerical results from Sec.|2|to investigate the deterioration of the endpoint lepton energy spectrum caused by 
beamstrahlung. The same analysis can be readily applied to smuon production in supcrsymmetry as well. 



2. ANALYTICAL APPROACHES TO BEAMSTRAHLUNG 

For physics purposes we need to know the energy distribution ijje{x,t) of the colliding electrons, where x = E/Eq 
is the electron energy as a fraction of the nominal beam energy Eq. Neglecting electrons from pair-production, the 
normalization condition on ipe is 

/ dxi;e{x,t) = 1. (1) 



*This talk was given by K. Kong, describing past and ongoing work performed in collaboration with the other authors. 
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Figure 1: (a) F{x,x') corresponds to the source and F{x",x) to the sink for the evolution of the spectral function ^e(x,t). 
(b) The final energy of an electron results from a convolution of beamstrahlung and bremsstrahlung (see Sec. 12.31 . 



Assuming that the emission of a photon takes place on an infinitesimally short time scale, the interference between 
successive radiation processes is negligible and the evolution of the electron spectral fmiction tp^ is described by the 
following rate equation (see Figure QJi): 

^^^^^^ = - j^^ dx"F{x\x)M^:-t)+ f^dx'F{x,x')M^\t) ■ (2) 

Here F(x,x') is the spectral function of radiation, i.e. F{x,x') is the transition probability (per unit time) for an 
electron at x' to move into the energy interval {x,x + dx). Obviously, F{x,x') = if a; > x'. Pulling out tjje{x,t), 
which is independent of x", outside the x" integral, leaves the integral quantity 

i/(x) = / dx" F{x" , x) , (3) 



Jo 

which represents the average number of photons radiated per unit time by the electron with an instantaneous energy 
X. The transition probability F{x,x') derived by Sokolov and Ternov P| is given by 

F{x,x') = '^m,v), (4) 
xx' 



1 



2^2 



n 



(5) 



where ^ = 3x'T/2, rj = k {1/x — 1/x') and the constant k = 2/(3T) is introduced for convenience. is the modified 
Bessel function and Vci is the number of photons per unit time calculated by the classical theory of radiation. By 
definition, this is also the limit of i'{x) for a; 0: 

5 

= v{x = 0) = — = T , (6) 

where Te is the classical electron radius, a is the fine structure constant, 70 = Eq/ {nieC?) and the dimcnsionlcss 
Lorentz invariant parameter 

6 a(Tz{crx + o-y) 

has a value specific to any given collider, as it depends on the total number of particles TV in a gaussian bunch and 
the rms sizes Gx, Cy, of the gaussian beam. 

2.1. Yokoya-Chen (YC) approximation 

It is clear that the Sokolov- Ternov (ST) spectral function Q is too complicated to allow a general analytic solution 
to Eq. (O. Therefore, one has to resort to some kind of an approximation. In 1989 Yokoya and Chen 0,13 suggested 
the following simple approximation to ((SJ: 

(8) 
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with the additional assumption 

r 

I'ix) = / dx" F{x, x") ~ Vci = constant. (9) 
Jo 

One can then solve the rate equation ^ by Laplace transformation and obtain an analytical solution; 



e 



(5(1 -x) + H7^l^^iy,,t) 

I — X 



(10) 



where r]x = k(1/x — 1) and 

^ i-X+ico °° 

h{u) = — exp{up~^/^ +p)dp = y] (11) 

27ri Ja-^oo n!r(n/3) 

with A > and < w < oo. The first term in Eq. Hl()|l represents the electron population unaffected by radiation. 
The n-th term in the Taylor expansion of the second term in 11U|I corresponds to the process of n-photon emission. 
The final electron energy distribution tpeix) due to beamstrahlung is obtained by suitable time averaging of the 
time-dependent solution (|10|) to account for the duration of bunch overlap, for example: 



2 /■'/■^ 1 

i^eix) = 7 / dt 1pe{x,t) = — 



(1 - e-^-')5{l -x) + ^^h{x) 

1 — X 



(12) 



where I = 2\/3az is the effective length of the incoming bunch and Nd = I'dl/'^ is the average number of photons 
radiated per particle during the entire collision of the e+e~ beams. The function h{x) in the second term is 

oo ri/3 

^(^)-E w( /oa (^ + l,^cO, (13) 
^-^ n'l [n/3) 

n— ^ ' ^ 

where 7(n, x) is the incomplete Gamma function 

7(n,a;)= C dye-yy''-\ (14) 
Jo 

The resulting solution 112|) is in good agreement with ABEL simulation data at least up to T ~ 0.44 

2.2. Consistent Yokoya-Chen (CYC) approximation 

Later on, M. Peskin pointed out 4| that the two assumptions and Q behind the YC approximation are 
incompatible with the normalization condition of the energy spectral function. To rectify the problem, Peskin 
proposed the so-called "consistent Yokoya-Chcn" approximation as an alternative to 

The rate equation Q can then be solved in exactly the same way and the resulting analytic solution is 

^e(x) = e-^- [six - 1) + ^-^^h{N,^l'^)^ , (16) 

where Nj ~ \f^ayVcii\ + T^/^)^^/^. This solution not only preserves the probability sum rule ^ but also agrees 
with Guinea Pig simulation for T = 0.104 (NLC500) and T = 0.299 (NLCIOOO) and numerically is not very different 
from the YC solution for T < 1 '2^. However, for large T, the two solutions arc significantly different, see Fig.|21 
and Sec. 13 below. 
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Figure 2: (a) Electron energy distributions due to beamstrahlung (dot-dashed, red), ISR (dashed, blue) and their convolution 
(solid), for an ILC with the parameters as shown. The YC and CYC approximations 1121 and 1161 are almost identical in 
this case, (b) Electron energy distributions due to beamstrahlung (dot-dashed) and ISR (dashed, green), for CLIC with the 
shown parameters. The YC (CYC) solution is shown in red (blue). 
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Figure 3: (a) The transition probability Fix, 1) as a function of x, using the exact expression (|KJ or the approximations JSJ 
and H15|l . for two different values of T. (b) The function vix) plotted for various values of T. 



2.3. Convolution with bremsstrahlung 

Another source of electron energy losses is bremsstrahlung, or initial state radiation (ISR), parametrized with its 
own energy distribution function ipisR The probability to end up with an electron of energy a; is a convolution 
of ISR and beamstrahlung (see Figs.^ and|21i): 

V'/SK * V'e = * V'/sfl = / dydz^pisR{y)ilJe{z)S{x-yz)= f ^ V/Sfl V'e(^) • (17) 



3. VALIDITY RANGE OF THE APPROXIMATE SOLUTIONS 

111 order to test tlie validity of the YC and CYC approximations, in Fig. we plot the transition probability 
F{x, 1) as a function of x, using the exact expression jSJ or the approximations (jS)) and (|15|) . for two different values 
of T. First of all, we notice that as T gets large, one enters the quantum regime of radiation, where emitting hard 
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Figure 4: (a) Muon energy spectrum from KK muon production in UED at an e^e~ collider with Ecm = 500 GeV and 
T = 0.107. The left (right) box-like distribution is due to SU(2)-singlet (SU(2)-doublet) /x^ production. The dotted lines 
ignore the beamstrahlung, with the effects from ISR shown in green. The red (blue) lines include beamstrahlung according to 
the YC (CYC) approximation, (b) The same as (a), but for R'^ = 500 GeV at CLIC with Ecm = 3000 GeV and T = 3.81. 



photons becomes more and more likely. The large tail at low x appearing in the case of CLIC could become an 
annoying spoiler for physics studies — an issue we will look at in more detail in the next Section. We also see that 
fixing the normalization problem of the YC approximation comes at a price: the improved CYC approximation is 
numerically further away from the exact ST spectral function Q and the discrepancy becomes more pronounced at 
large T. Conversely, while the YC approximation docs a better job in approximating it not only suffers from the 
normalization problem, but also relies on the approximation v{x) « Vd = const, which fails at large T, as evidenced 
in Fig.|2|3. There we plotted v{x) and the corresponding Vd for different values of T. We see that for small T the 
v{x) function (j2J) is pretty constant and the YC approximation is valid, but for the larger T values which arc relevant 
at CLIC, is far from being a constant, and the YC approximation will fail. 

Fig. 13 reveals that in order to get a reliable parametrization of beamstrahlung for large T (i.e. at very high energy 
e+e~ colliders such as CLIC), we need to improve on the two existing approximations. We tried, but were unable 
to find a better approximation to Q, which would still allow for an analytic solution in closed form. We therefore 
resorted to a numerical solution of the original rate equation |5J) with the full ST formula This solution was then 
used for the physics study presented in Scc.0| 



4. MUON ENERGY SPECTRUM FROM KK-MUON PRODUCTION 

In order to study the effect of beamstrahlung on physical observables. we consider the production of level 1 Kaluza- 
Klein muons /if in UED 0,0], which then decay to the lightest KK particle (LKP), in this case the KK partner 71 of 
the photon. The signature is two opposite sign muons and missing energy, which is analogous to smuon production 
in supersymmetry models with stable neutralino LSP. At hadron colliders, the two scenarios can be confused 
but lepton colliders allow for an easy and straightforward discrimination . 

Two important observables arc the lower. Emin, and upper, Emax, endpoints of the muon energy spectrum (see 
Fig. 2J). They can be used to extract the masses of the particles involved in the decay as follows: 

1 / A/^ \ 

Eraa./^^n^l^M^Al-J^U{l±^3), (18) 

where M^^ and M^^ are the KK muon and LKP masses, respectively. 7 = 1/(1 — and /3 is the fii boost. 

Fig. confirms that at low Ecm the effects from beamstrahlung are rather small, and furthermore, the YC and 
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Figure 5: (a) Electron energy distributions resulting from the numerical solution (green) or the YC (red) and CYC (blue) 
approximations, for T = 3.81 and T = 0.107. (b) The same as Fig.^J), but including the exact solution for beamstrahlung. 



CYC approximations are in agreement. However, Fig. 2Jd reveals that at higher energy colliders, beamstrahlung 
tends to wash out the endpoint, the degree of deterioration being quite sensitive to the approximation used. 

Since we already expect neither approximations to be reliable at large Ecm (where T is also large), we solved 
Eq. |5J) numerically. The result, together with the YC and CYC approximations, is shown in Fig. Et, for T = 
0.107 (TESLA500) and T = 3.81 (CLIC3000) Notice how the improved CYC approximation considerably 

overestimates the effect of beamstrahlung at low x. We then implemented the numerical solution in the CompHEP 
event generator [ll] |. The resulting true muon energy distribution including both types of beam energy losses is 
shown in Fig. Oj. We see that beamstrahlung still degrades the kinematic endpoint, although not as much as the 
CYC and YC approximations would suggest. 
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